High number of SARS-CoV-2 persistent infections uncovered through genetic analysis of samples from a large community-based surveillance study (preprint)

Persistent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections may act as viral reservoirs that could seed future outbreaks 1-5, give rise to highly divergent lineages 6-8, and contribute to cases with post-acute Coronavirus disease 2019 (COVID-19) sequelae (Long Covid) 9,10. However, the population prevalence of persistent infections, their viral load kinetics, and evolutionary dynamics over the course of infections remain largely unknown. We identified 381 infections lasting at least 30 days, of which 54 lasted at least 60 days. These persistently infected individuals had more than 50% higher odds of self-reporting Long Covid compared to the infected controls, and we estimate that 0.09-0.5% of SARS-CoV-2 infections can become persistent and last for at least 60 days. In nearly 70% of the persistent infections we identified, there were long periods during which there were no consensus changes in virus sequences, consistent with prolonged presence of non-replicating virus. Our findings also suggest reinfections with the same major lineage are rare and that many persistent infections are characterised by relapsing viral load dynamics. Furthermore, we found a strong signal for positive selection during persistent infections, with multiple amino acid substitutions in the Spike and ORF1ab genes emerging independently in different individuals, including mutations that are lineage-defining for SARS-CoV-2 variants, at target sites for several monoclonal antibodies, and commonly found in immunocompromised patients 11-14. This work has significant implications for understanding and characterising SARS-CoV-2 infection, epidemiology, and evolution.

Viral burdens are associated with age and viral variant in a population-representative study of SARS-CoV-2 that accounts for time-since-infection related sampling bias. (preprint)

In this study, we evaluated the impact of viral variant, in addition to other variables, on within-host viral burdens, by analysing cycle threshold (Ct) values derived from nose and throat swabs, collected as part of the UK COVID-19 Infection Survey. Because viral burden distributions determined from community survey data can be biased due to the impact of variant epidemiology on the time-since-infection of samples, we developed a method to explicitly adjust observed Ct value distributions to account for the expected bias. Analysing the adjusted Ct values using partial least squares regression, we found that among unvaccinated individuals with no known prior infection, the average Ct value was 0.94 lower among Alpha variant infections, compared those with the predecessor strain, B.1.177. However, among vaccinated individuals, it was 0.34 lower among Delta variant infections, compared to those with the Alpha variant. In addition, the average Ct value decreased by 0.20 for every 10 year age increment of the infected individual. In summary, within-host viral burdens are associated with age, in addition to the interplay of vaccination status and viral variant.

Lineage replacement and evolution captured by the United Kingdom Covid Infection Survey (preprint)

The Office for National Statistics COVID-19 Infection Survey is a large household-based surveillance study based in the United Kingdom. Here, we report on the epidemiological and evolutionary dynamics of SARS-CoV-2 determined by analysing sequenced samples collected up until 13th November 2021. We observed four distinct sweeps or partial-sweeps, by lineages B.1.177, B.1.1.7/Alpha, B.1.617.2/Delta, and finally AY.4.2, a sublineage of B.1.617.2, with each sweeping lineage having a distinct growth advantage compared to their predecessors. Evolution was characterised by steady rates of evolution and increasing diversity within lineages, but with step increases in divergence associated with each sweeping major lineage, leading to a faster overall rate of evolution and fluctuating levels of diversity. These observations highlight the value of viral sequencing integrated into community surveillance studies to monitor the viral epidemiology and evolution of SARS-CoV-2, and potentially other pathogens, particularly as routine PCR testing is phased out or in settings where large-scale sequencing is not feasible.

Highly-Sensitive Lineage Discrimination of SARS-CoV-2 Variants through Allele-Specific Probe Polymerase Chain Reaction (preprint)

Introduction: Tools to detect SARS-Coronavirus-2 variants of concern and track the ongoing evolution of the virus are necessary to support ongoing public health efforts and the design and evaluation of novel COVID-19 therapeutics and vaccines. Although next-generation sequencing (NGS) has been adopted as the gold standard method for discriminating SARS-CoV-2 lineages, alternative methods may be required when processing samples with low viral loads or low RNA quality. Methods: An allele-specific probe polymerase chain reaction (ASP-PCR) targeting lineage-specific single nucleotide polymorphisms (SNPs) was developed and used to screen 1,082 samples from two clinical trials in the United Kingdom and Brazil. Probit regression models were developed to compare ASP-PCR performance against 1,771 NGS results for the same cohorts. Results: Individual SNPs were shown to readily identify specific variants of concern. ASP-PCR was shown to discriminate SARS-CoV-2 lineages with a higher likelihood than NGS over a wide range of viral loads. Comparative advantage for ASP-PCR over NGS was most pronounced in samples with Ct values between 26-30 and in samples that showed evidence of degradation. Results for samples screened by ASP-PCR and NGS showed 99% concordant results. Discussion: ASP-PCR is well-suited to augment but not replace NGS. The method can differentiate SARS-COV-2 lineages with high accuracy and would be best deployed to screen samples with lower viral loads or that may suffer from degradation. Future work should investigate further destabilization from primer:target base mismatch through altered oligonucleotide chemistry or chemical additives.

Efficacy of ChAdOx1 nCoV-19 (AZD1222) vaccine against SARS-CoV-2 lineages circulating in Brazil; an exploratory analysis of a randomised controlled trial (preprint)

Background Emerging evidence shows the substantial real-world impact of authorised vaccines against COVID-19 and provides insight into the potential role of vaccines in curbing the pandemic. However, there remains uncertainty about the efficacy of vaccines against different variants of the virus. Here we assessed efficacy of ChAdOx1 nCoV-19 (AZD1222) against lineages of SARS-CoV-2 circulating in Brazil from June 2020 until early 2021. Methods Participants aged 18 and above were enrolled into a randomised phase 3 trial of ChAdOx1 nCoV-19 vaccine against symptomatic SARS-CoV-2 infection. Participants received two doses of ChAdOx1 nCoV-19 or control (1st dose: Men ACWY vaccine, 2nd dose: normal saline). Nasopharyngeal and oropharyngeal swabbing was performed if participants developed symptoms of COVID-19 (cough, shortness of breath, fever >37.8°C, ageusia, anosmia). Swabs were tested by nucleic acid amplification (NAAT) for SARS-CoV-2, sequenced, and viral load determined. For those samples where a genotype could not be ascertained from sequencing, allele specific PCR was performed. The efficacy analysis included symptomatic COVID-19 in seronegative participants with a NAAT positive swab more than 14 days after a second dose of vaccine. Participants were unblinded after the vaccine was authorised for use, and the control participants offered vaccination. Infections occurring after unblinding were excluded from analysis. Vaccine efficacy was calculated as 100% x (1 – relative risk (RR)), where RR was estimated from a robust Poisson model. The trial is registered at ISRCTN89951424. Findings 9433 participants were eligible for inclusion in the pre-specified primary efficacy population, having reached more than 14 days after a second dose of ChAdOx1 nCoV-19, of whom 307 were NAAT+, in this post-hoc analysis. From June 2020 to February 2021, the two most frequently identified lineages were P.2 (N=153) and B.1.1.28 (N=49). P.1 emerged during the study (N=18) but became dominant only after study unblinding. Viral loads were highest amongst those with P.1 infection. Vaccine efficacy (VE) for B.1.1.33 (88.2%, 95%CI 5, 99), B.1.1.28 (73%, 95% CI, 46, 86), P.2 (69% 95% CI, 55, 78) and P.1 (64%, 95% CI, -2, 87) was estimated. In participants who had received two doses of vaccine, one COVID-19 hospitalisation occurred in the ChAdOx1 nCoV-19 group and 18 in the control group, with VE against hospitalisation 95% (95% CI 61, 99). There were 2 COVID-19 deaths in the control group and none in the vaccine group. Interpretation ChAdOx1 nCoV-19 provides high efficacy against hospitalisation, severe disease and death from COVID-19 in Brazil and there is strong evidence of protection being maintained against P.2, despite the presence of the spike protein mutation E484K. Real world effectiveness studies are ongoing in Brazil to further establish protection against P.1 and other emerging variants.

A blood atlas of COVID-19 defines hallmarks of disease severity and specificity (preprint)

Treatment of severe COVID-19 is currently limited by clinical heterogeneity and incomplete understanding of potentially druggable immune mediators of disease. To advance this, we present a comprehensive multi-omic blood atlas in patients with varying COVID-19 severity and compare with influenza, sepsis and healthy volunteers. We identify immune signatures and correlates of host response. Hallmarks of disease severity revealed cells, their inflammatory mediators and networks as potential therapeutic targets, including progenitor cells and specific myeloid and lymphocyte subsets, features of the immune repertoire, acute phase response, metabolism and coagulation. Persisting immune activation involving AP-1/p38MAPK was a specific feature of COVID-19. The plasma proteome enabled sub-phenotyping into patient clusters, predictive of severity and outcome. Tensor and matrix decomposition of the overall dataset revealed feature groupings linked with disease severity and specificity. Our systems-based integrative approach and blood atlas will inform future drug development, clinical trial design and personalised medicine approaches for COVID-19.

Virological and serological characterization of critically ill patient with COVID-19 in the UK: a special focus on variant detection (preprint)

Background. Treatment of COVID-19 patients with convalescent plasma containing neutralising antibody to SARS-CoV-2 is under investigation as a means of reducing viral loads, ameliorating disease outcomes, and reducing mortality. However, its efficacy might be reduced in those infected with the emerging B.1.1.7 SARS-CoV-2 variant. Here, we report the diverse virological characteristics of UK patients enrolled in the Immunoglobulin Domain of the REMAP-CAP randomised controlled trial. Methods. SARS-CoV-2 viral RNA was detected and quantified by real-time PCR in nasopharyngeal swabs obtained from study subjects within 48 hours of admission to intensive care unit. Antibody status was determined by spike-protein ELISA. B.1.1.7 strain was differentiated from other SARS-CoV-2 strains by two novel typing methods detecting the B.1.1.7-associated D1118H mutation with allele-specific probes and by restriction site polymorphism (SfcI). Findings. Of 1260 subjects, 90% were PCR-positive with viral loads in nasopharyngeal swabs ranging from 72 international units [IUs]/ml to 1.7x10^11 IU/ml. Median viral loads were 45-fold higher in those who were seronegative for IgG antibodies (n=314; 28%) compared to seropositives (n=804; 72%), reflecting in part the latter group's possible later disease stage on enrolment. Frequencies of B.1.1.7 infection increased from early November (<1%) to December 2020 (>60%). Anti-SARS-CoV-2 seronegative individuals infected with wild-type SARS-CoV-2 had significantly higher viral loads than seropositives (medians of 1.2x10^6 and 3.4 x10^4 IU/ml respectively; p=2x10^-9). However, viral load distributions were elevated in both seropositive and seronegative subjects infected with B.1.1.7 (13.4x10^6 and 7.6x10^6 IU/ml; p=0.18). Interpretation. High viral loads in seropositive B.1.1.7-infected subjects are consistent with increased replication capacity and/or less effective clearance by innate or adaptive immune response of B.1.1.7 strain than wild-type. As viral genotype was associated with diverse virological and immunological phenotypes, metrics of viral load, antibody status and infecting strain should be used to define subgroups for analysis of treatment efficacy.

Early Analysis of a potential link between viral load and the N501Y mutation in the SARS-COV-2 spike protein (preprint)

A new variant of SARS-CoV-2 has emerged which is increasing in frequency, primarily in the South East of England (lineage B.1.1.7 (1); VUI-202012/01). One potential hypothesis is that infection with the new variant results in higher viral loads, which in turn may make the virus more transmissible. We found higher (sequence derived) viral loads in samples from individuals infected with the new variant with median inferred viral loads were three-fold higher in individuals with the new variant. Most of the new variants were sampled in Kent and Greater London. We observed higher viral loads in Kent compared to Greater London for both the new variant and other circulating lineages. Outside Greater London, the variant has higher viral loads, whereas within Greater London, the new variant does not have significantly higher viral loads compared to other circulating lineages. Higher variant viral loads outside Greater London could be due to demographic effects, such as a faster variant growth rate compared to other lineages or concentration in particular age-groups. However, our analysis does not exclude a causal link between infection with the new variant and higher viral loads. This is a preliminary analysis and further work is needed to investigate any potential causal link between infection with this new variant and higher viral loads, and whether this results in higher transmissibility, severity of infection, or affects relative rates of symptomatic and asymptomatic infection Document Description and PurposeThis is an updated report submitted to NERVTAG in December 2020 as part of urgent investigations into the new variant of SARS-COV-2 (VUI-202012/01). It makes full use of (and is restricted to) all sequence data and associated metadata available to us at the time this original report was submitted and remains provisional. Under normal circumstances more genomes and metadata would be obtained and included before making this report public. We will update this preprint when more genomes and metadata are available and before submitting for peer review.

SARS-CoV-2 genomics beyond the consensus sequence: evidence for circulating mixed viral populations (preprint)

We gratefully acknowledge the UK COVID-19 Genomics Consortium (COG UK) for funding, and Public Health Wales / Cardiff University and MRC-University of Glasgow Centre for Virus Research for making their COG-UK sequence data publicly available. COG-UK is supported by funding from the Medical Research Council (MRC) part of UK Research & Innovation (UKRI), the National Institute of Health Research (NIHR) and Genome Research Limited, operating as the Wellcome Sanger Institute. The research was supported by the Wellcome Trust Core Award Grant Number 203141/Z/16/Z with funding from the NIHR Oxford BRC. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health. We are deeply grateful to Robert Esnouf and the BMRC Research Computing team for unfailing assistance with computational infrastructure. We also thank Benjamin Carpenter and James Docker for assistance in the laboratory, and Lorne Lonie, Maria Lopopolo, Chris Allen, John Broxholme and the WHG high-throughput genomics team for sequencing and quality control. The HIV clone p92BR025.8 was obtained through the Centre For AIDS Reagents from Drs Beatrice Hahn and Feng Gao, and the UNAIDS Virus Network (courtesy of the NIH AIDS Research and Reference Reagent Program). KAL is supported by The Wellcome Trust and The Royal Society (107652/Z/15/Z). MH, LF, MdC, GMC, NO, LAD, DB, CF and TG are supported by Li Ka Shing Foundation funding awarded to CF. PS is supported by a Wellcome Investigator Award (WT103767MA). SummarySARS-CoV-2, the causative agent of COVID-19, emerged in late 2019 causing a global pandemic, with the United Kingdom (UK) one of the hardest hit countries. Rapid sequencing and publication of consensus genomes have enabled phylogenetic analysis of the virus, demonstrating SARS-CoV-2 evolves relatively slowly1, but with multiple sites in the genome that appear inconsistent with the overall consensus phylogeny2. To understand these discrepancies, we used veSEQ3, a targeted RNA-seq approach, to quantify minor allele frequencies in 413 clinical samples from two UK locations. We show that SARS-CoV-2 infections are characterised by extensive within-host diversity, which is frequently shared among infected individuals with patterns consistent with geographical structure. These results were reproducible in data from two other sequencing locations in the UK, where we find evidence of mixed infection by major circulating lineages with patterns that cannot readily be explained by artefacts in the data. We conclude that SARS-CoV-2 diversity is transmissible, and propose that geographic patterns are generated by transient co-circulation of distinct viral populations. Co-transmission of mixed populations could open opportunities for resolving clusters of transmission and understanding pathogenesis.